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Abstract: The problem of Bose-Einstein condensation for a relativistic ideal gas on 
a 3+1 dimensional manifold with a hyperbolic spatial part is analyzed in some 
detail. The critical temperature is evaluated and its dependence of curvature 

R -1 ' is pointed out. 
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1 Introduction 

Bose-Einstein condensation for a non relativistic ideal gas has a long history |I|. The 
physical phenomenon is well described in many text books (see for example ref . || ) 
and a rigorous mathematical discussion of it was given by many authors]^, [|. The 
generalization to a relativistic idel Bose gas is non trivial and only recently has been 
discussed in a series of papers^], ^|, |7[. 

*Email: cognola@itnvax.cineca.it, vanzo@itncisca 



It is well known that in the thermodynamic limit (infinite volume and fixed 
density) there is a phase transition of the first kind in correspondence of the critical 
temperature at which the condensation manifests itself. At that temperature, the 
first derivative of some continuous thermodynamic quantities has a jump. If the 
volume is keeped finite there is no phase transition, nevertheless the phenomenon 
of condensation still occurs, but the critical temperature in this case is not well 
defined. 

For manifolds with compact hyperbolic spatial part of the kind H /T, V be- 
ing a discrete group of isometries for the N-dimensional Lobachevsky space H N , 
zero temperature effects as well as finite temperature effects induced by non-trivial 
topology, have been recently studied in some detail || |, [10], [II], [12], ?, [H], [15], [16]. 



To our knowledge, a similar analysis has not yet been carried out for non compact 
hyperbolic manifolds. 

Hyperbolic spaces have remarkable properties. For example, the continuous spec- 
trum of the Laplace-Beltrami operator has a gap determined by the curvature ra- 
dius of H N , implying that massless fields have correlation functions exponentially 
decreasing at infinity (such a gap is not present for the Dirac operator). For that 
reason, H 4 was recently proposed as an excellent infrared regulator for massless 
quantum field theory and QCD|L7|]. Critical behaviour is even more striking. In two 
flat dimension vortex configurations of a complex scalar field, the XY model for He 4 
films have energy logarithmically divergent with distance, while on H 2 it is finite. 
This implies that the XY model is disordered at any finite temperature on H 2 . Even 
quantum mechanics on H 2 has been the subject of extensive investigations [15] . The 
manifold H 4 is also of interest as it is the Euclidean section of anti-de Sitter space 
which emerges as the ground state of extended supergravity theories. The stress 
tensor on this manifold has been recently computed for both boson and fermion 
fields using zeta-function methods JIU[. 

In the present paper we shall discuss finite temperature effects and in particular 
the Bose-Einstein condensation for a relativistic ideal gas in a 3+1 dimensional 
ultrastatic space-time M. = RxH 3 . We focus our attention just on H 3 , because such 
a manifold could be really relevant for cosmological and astrophysical applications. 



To this aim we shall derive the thermodynamic potential for a charged scalar field of 
mass m on Ai, using zeta function, which on H 3 is exactly known. We shall see that 
the thermodynamic potential has two branch points when the chemical potential \i 
riches ±u; , ou 2 — K + m 2 being the lower bound of the spectrum of the operator 
L m = —A + m 2 and — k the negative constant curvature of H 3 . The values ±uj 
will be riched by \x — /u(T) of course for T = 0, but also for T = T c > 0. This is the 
critical temperature at which the Bose gas condensates. 

The paper is organized as follows. In section 2 we study the elementary properties 
of the Laplace-Beltrami operator on H 3 ; in particular we derive its spectrum and 
build up from it the related zeta-function. In section 3 we briefly recall how zeta- 
function can be used in order to regularize the partition function and we derive 
the regularized expression for the thermodynamic potential. In section 4 we discuss 
the Bose-Einstein condensation and derive the critical temperatures in both the 
cases of low and high temperatures. In section 5 we consider in detail the low and 
high temperature limits and derive the jump of the first derivative of the specific 
heat. The paper end with some considerations on the results obtained and some 
suggestions for further developments. 

2 The spectrum and the zeta function of Laplace- 
Beltrami operator on H 3 

For the aims of the present paper, the 3-dimensional Lobachevsky space H 3 can be 
seen as a Riemannian manifold of constant negative curvature —k, with hyperbolic 
metric dl 2 = dg 2 + sinh g{d"& 2 -\-s\v? •& d"d 2 ) and measure dVt = sinh gdgdT,, oE being 
the measure on S 2 . For convenience, here we normalize the curvature — k to — 1. In 
these coordinates, the Laplace-Beltrami operator A reads 

d 2 d 1 

A = — + 2coth£— + — r^A 52 (1) 

og 2 og sinh g 

A52 being the Laplacian on the sphere. 

The zeta function related to —A, can be easily computed if one knows the density 
of states. For our aim, only states with zero angular momentum will be important, 



so, in order to derive it, it is sufficient to study radial wave functions of —A, that 

is solution of equation 

d 2 u , du .„. 

— -= + 2 coth — + Xu = (2) 

dg 2 dg 

which reduces to 

with the substitution u — vj sinh g. The bounded solution of the latter equation is 
trivially given by v = A sin ug, v being related to the eigenvalue by A = v 2 + 1. Such 
a solution remains bounded at infinity provided v stays real, or A > 1. Thus the 
continuum spectrum of —A on H 3 has a lower bound at A = 1 (or A = k in standard 
unities). Note that in general this is not true on the compact manifold H 3 /T, which 



has been considered in ref.[pl|. We also notice that the wave operator propagates 

the field excitations on the light cone, hence the gap should not be interpreted as a 

physical mass. 

For the radial solution normalized to 1 at the origin one gets 

/ x smug 
uM = —r- (4) 

v sinh q 

Now, the L 2 (dQ) scalar product for u u (g) is 



(u u ,u u i) — / smv g sin v' gdg = — —5{v — v') (5) 

vv' h v l 

from which the density of states giv) = Vv 2 /2it 2 directly follows. As usual, we have 

introduced the large, finite volume V to avoid divergences. When possible, the limit 

V — > oo shall be understood. 

At this point the computation of zeta function is straightforward. As we shall see 

in the following, what we are really interested in, is the zeta function related to the 

operators Q± = L^ 2 ± /j,. The eigenvalues of L m are <jj 2 {y) = v 2 + a 2 = v 2 + k + m 2 , 

then we get 



C(s;Q±) = — 2 j o {u{v)±n)- s v 2 dv 

• S|(S ;J))V -F( S + M- 3; *-!;**) (6) 



(47r) 3 /2r( s -l/2)(2a 



,s-3 V" i "->" ")" 2' 2a 



where F(a, /3; 7; 2;) is the hypergeometric function. For its properties and its integral 



representations see for example ref.[20[. It has to be noted that eq.(|J) is the very 
same one has on a flat space for a massive field with mass equal to a. Here in fact, 
the curvature plays the role of an effective mass. 

As we see from eq.fl^), the zeta function related to the pseudo-differential oper- 
ators Q± has simple poles at the points s n = 3, 2, 1, — 1, — 2, —3, . . . with residues 
b n (±jj) = Res (C(s; Q±), s n ) given by 



V Va V 

— ; &2 (±^) = T _ ; 6 1 (± A *) = — 



& 3 ( ±/i ) = _; U±») = T-£; h{±v) =-- —{2ii 2 - a 2 )- (7) 



6_ n (±/i) = c. n F(-n + 1, -n- 3; -n - 1/2;^); (8) 

(-l) n nV r (2a) n+3 
°~ n ~ (47r) 3 / 2 r(-n-l/2)r(n + 4) ; n - ; 

Moreover, using the following property of the hypergeometric function: 

i-F(a, /9; 7; «) = — F(a + 1, /? + 1; 7 + M (9) 

az 7 

one easily sees that all fc_ n (±/x) are odd or even polinomials of /i, according to 
whether n is even or odd; this means that 6_ n (/i) = (— l) n+1 6_ n (— //) and the function 
C(Q+) + C(Q-) nas simple poles at the points s = —2k + 3 with residues given by 
26_ 2 fc+3(/i) and simple zeros for s = —2k (k = 0, 1,2, . . .). For the special values 
ji = ±a we immediately get 6_„(a) = (— l) n+1 6_„(— a) = c_ n . Finally from eq.@ 
we also derive the recurrence relations 

d^b_ n (±fi) = d=nfo_ (n _i) (10) 

3 The thermodynamic potential 

In order to define the grand canonical partition function Z(j3, //) = e - ^ n ^'^, f2(/3, /x) 
being the thermodynamic potential, we suppose the charged scalar field to be in 
thermal equilibrium at finite temperature T = 1/(3. Then, according to ref. [21], one 
has 

Z(J3,fi)= I V0V(f)e- ti ' dT Jt^VadV ^ 

J<t>(0,x)=<t>(l3,x) 



where L M = — (d T — /i) 2 + L m and the Wick rotation r = ix° has to be understood. 
In eq.flTTD the integration has to be taken over all fields <f>(r, x a ) with /3-periodicity 
with respect to r. 

The eigenvalues of the whole operator L M , say \i n ^ read 

Hn,» = [ — + ifj.\ +u;(v) 2 n = 0,±l,±2,... (12) 

Regularizing the partition function in eq.([TT|) in terms of zeta-function ((s; L M ) 
related to the operator L M [Q one easily gets 



n(P,ft) = ilogdet(r 2 L M ) = -iC'(0;r 2 L M ) 

= -i[log^C(0;i M ) + C'(0;^)] (13) 

£ being an arbitrary normalization parameter coming from the scalar path-integral 
measure. Note that £, which has the dimensions of a mass, is necessary in order 
to keep the zeta-function dimensionless for all s. The finite temperature and yu 
dependent part of the thermodynamic potential does not suffer of the presence of 
such an arbitrary parameter. On the contrary, £ enters in the regularized expression 
of vacuum energy and this creates an ambiguity [12], which is proportional to the 



heat kernel expansion coefficient K N (L m ) related to L m (in general, K N (L m ) ^ 0). 
When the theory has a natural scale parameter, like the mass of the particle or the 
constant curvature of the manifold, the ambiguity can be removed by an " ad hoc" 



choice of £ 23 



Here we would like to study the behaviour of thermodynamic quantities, then we 
are only interested in the /i and T dependent part of the thermodynamic potential; 
that is a well defined quantity, which does not need regularization. To compute it, 
it is not necessary to use all the analytic properties of zeta function (for a careful 



derivation of vacuum energy see for example refs. |]I2|, |l4f). Then we can proceed 
in a formal way and directly compute log det L M disregarding the vacuum energy 
divergent term. 

First of all we observe that 



JT log(u 2 + (2™/ (3 + ifi) 2 ) 



diO 7 



uj 2 + (2im//3 + ifif 



E 

?1= — oo 

f j- [coth I (w + fi) + coth | {uj - //)] dw^U) 
- log [l - e-^+^l - log [l - e-^-^] - (3u 



Using eq . (|T3|) , recalling that u 2 = v 2 + a 2 and by integrating over v with the state 
density that we have derived in the previous section, we get the standard result 



ftOM 



--Tr logL M 



E(P,im)+E v 



V 



2n 2 /3 
V 

+ 2^ 



lot 



-P(w(v)+li) 



+ log 



-P{u{v)-y) 



}u 2 du (15) 



ijj{y)v 2 dv 



The last term is divergent and represents the non regularized vacuum energy that we 
throw away. The rest of the equation is exactly the \x and T dependent part of the 
thermodynamic potential from which we can derive all thermodynamic quantities. 
0(/3, fi) as a function of the complex parameter fi, has branch points when fi 2 is equal 
to some oj 2 (v). Thus the physical values of fi must be restricted by \fi\ < u = a. 
As we shall see in the next section, the value \fi\ — a can be riched for a critical 
temperature at which Bose-Einstein condensation takes place. 

In order to get low and high temperature expansions, it is convenient to get 
different representations of E(/3, fi), which can be derived in a rigorous way starting 
from eq. (|T3|) . For a compact manifold H 3 /T they have been derived in ref. |15[, to 
which we refer the reader for more details. Here we just write down the final results, 
observing that the expressions one has on H 3 are the same one has on H 3 /T when 
topological contributions are dropped out. Then we have 



E(/3,fi) 
E{(3,fi) 
E{(3,fj) 



aV» K 2 (an(3) 
— — 2_^ cosh np/j, 



TT^ 



2R2 



n=l 



n 2 (3 



(16) 



717 



;E 



/' 



2d 



' T{2n + 1) J Re 



[ r( S +2n-l)C R (s)C(^f^; L m ) f3~ s ds(17) 

J tie s=c 



- f T(s)Us + 1) [C(s; Q+) + C(a; Q-)\ P~ (s+1) da (18) 

ni J tie s=c 



where K2 is the modified Bessel function, c is a sufficiently large real number and 
(r(s) is the usual Riemann zeta-function. The integral representations (|TTD and 
(|l8l), which are valid for \p\ < a, are useful for high temperature expansion. On the 
contrary, the representation (|T6|) in terms of modified Bessel functions is more useful 
for the low temperature expansion, since the asymptotics of K v is well known. 

4 Bose-Einstein condensation 

In order to discuss Bose-Einstein condensation we have to analyze the behaviour of 
the charge density 

dQ(V,(3,z) _ 
P = z ^ = f{z)-f(l/z) (19) 

in the infinite volume limit. Here 

/(*) = Ev? — I v ( 2 °) 

y 1 V(exp fiujj - z) 

and the activity z = expftp has been introduced. The Uj in the sum are meant 
to be the Dirichlet eigenvalues for any normal domain V C H 3 . That is, V is a 
smooth connected submanifold of H 3 with non empty piecewise C°° boundary. By 
the infinite volume limit we shall mean that a nested sequence of normal domains 
Vk has been choosen together with Dirichlet boundary conditions and such that 
U Vk = H 3 . The reason for this choice is the following theorem due to Mac Kean 
(see for example @): 

- if u k denotes the smallest Dirichlet eigenvalue for any sequence of normal 
domains Vk filling all of H 3 then u k > a and lim^oo uJok = 0,. (Although the 
above inequality is also true for Neumann boundary conditions, the existence 
of the limit in not assured to the authors knowledge). 

Now we can show the convergence of the finite volume activity z^ to a limit point 
z as k — > 00. To fix ideas, let us suppose p > 0: then z k e (l,exp/3u; fc). Since 
p(V, /3, z) is an increasing function of z such that p(V, /3, 1) = and p(V, j3, 00) = 00, 
for each fixed Vk there is a unique Zk(p, 0) € (1, j3exp u k) such that p = p{Vk, /3, zu). 



By compactness, the sequence Zk must have at least one fixed point z and as u k -^ ° 2 
as k goes to infinity, by Mc Kean theorem, z G [1, exp (3a\. 

From this point on, the mathematical analysis of the infinite volume limit exactly 
parallels the one in flat space for non relativistic systems, as it is done in various 
references p5| , |26| , |J. In particular, there is a critical temperature T c over which 
there are no particles in the ground state. T c is the unique solution of the equation 

sinh (3a f°° v 2 dv 

[21) 



2vr 2 Jo cosh (3a Jl + (u/a) 2 - cosh (3a 



For T > T c we have |/i| < a, the charge density of the ground state g — and 
p is given by 



Q = 



Vd/j, Jo 



1 



& 



P(u>(v)+n) _ I e p{w{y)-n) 



-av 



2tt 2 



sinh/3/x f°° v l dv 



f22l 



2tt 2 Jo cosh (3uj(v) — cosh (3 pi 

For T < T c , \fi\ remains equal to a, the charge density g ex of the particles in the 
excited states is given by eq.(^T|) and the charge density of the particles in the 
ground state is non vanishing and given by g — g — g ex . That is, below T c we have 
Bose-Einstein condensation. 

So far we only have determined \x in the condensation region where in fact it 
is equal to a. In the uncondensed region it is determined by eq.(^) as an implicit 
function of g and T. For massive particles, it looks like the one, one has for the flat 
casef|[], with the only difference that the mass m is replaced by a = y/n + m 2 . The 
very difference between flat and hyperbolic spaces occurs for massless particles. We 
shall return on this important point in a moment. 

Solutions of eq. (^T|) can be easily obtained in the two cases (3a 3> 1 and (3a <C 1 
(in the case of massive bosons these correspond to non relativistic and ultrarela- 
tivistic limits respectively). We have in fact 

T 3 r°° Tx 2 dx 



2vr 2 Jo e x2 / 2a - 1 
aT 2 r°° x 2 dx 



27r 2 Jo coshx 



'oT\ 3 / 2 


(3a > 1 


(23) 


aT 2 


(3a < 1 


(24) 



T - 


2vr ( g 


J-c 


a Ur(3/2) 


T - 

J- c 


■(*)" 



from which we get the corresponding critical temperatures 

2/3 

| ; (3a > 1 (25) 

(5a < 1 (26) 

and the corresponding charge densities of particles in the ground state 

Qo = g[l- {T/T c f' 2 ] ; pa > 1 (27) 

q = q[1- [T/T c f] ■ (3a < 1 (28) 

In the flat case, that is for k = 0, eqs. (|25|) and (p7|) reduce to the results given in 
many textbooks (see for example ref. 0), while eqs.(^) and ( p8|) agree with the ones 
given in ref.||. 

It has to bo noted that for massless bosons, the condition \y\ < a does not 
require // = like in the flat space, because a > also for massless particles. This 
implies that the critical temperature is always finite and so, unlike in the flat case, q 
is always different from g . As has been noticed in ref.|J, on a flat manifold, the net 
charge of massless bosons resides in the Bose-Einstein condensed ground state. This 
never happens if the spatial manifold is H 3 . In fact, because of curvature, massive 
and massless bosons have a similar behaviour. 

As it is well known, at the critical temperature, continuous thermodynamic quan- 
tities may have a discontinuous derivative (first order phase transition). This derive 
from the fact that the second derivative of the chemical potential \x is discontinues 
for T = T c . Of course, // = dfi/dT = for T < T c , since in this region \x is con- 
stant. By considering the constant charge density pasa function of temperature 
and taking the total derivative with respect to T of eq.(^), we see that 

"' = -¥¥\ < 29 > 

and since d^g diverges for /x = a we obtain /x'(T c + ) = 0. This is not the case of 
//'. In fact we shall see that //'(T c + ) is different from zero and therefore n"(T) is a 



discontinuous function of temperature. This implay that the first derivative of the 
specific heat Cy has a jump for T = T c given by 



dC- 



v 



dT 



da 



v 



dT 






- t 6 V(t; 



dg(T, fj) 



T=T? 



dT 



(30) 



T=Ti 



U(T, fi) being the internal energy, which can be derived by means of equation 



(31) 



To go further in the computation, it is convenient to integrate eq.(17) in order 
to obtain a series representation of E(/3, //) valid for T > T c . This is the topic of the 
next section. 



5 Low and high temperature expansions 

The low temperature expansion can be easily obtained by using eq. (|i~6|) and recalling 
that for real values of s the asymptotics for the modified Bessel function reads 

r(fc + 5/2) 



Ko(s) 



£«-E 



2s ^ T(k + l)T(-k + 5/2) 
Then, for small T we have 



(2 S )- 



(32) 



£(/?,/!) 



a 4 V 



E 



-nf){a—\ii\) 



5/2 



and deriving with respect to \x 

e(T,/*)^-( — ) £ 



an(3 



E 

fc=0 



r(fc + 5/2) 



V{k + l)T(-k + 5/2) 



(2an/3)~ k (33) 



V27T/ 



n=l 



/r 



3/2 



(34) 



where only the leading terms have been taken into account. Setting |/i| = a we again 
obtain the charge density (|23|). By deriving eq.(|34"D we obtain 



aA?c 



2C 2 
Tl2 : 



™<9£ 3 

T 9T = 2^ 



(35) 



where A = 2.363 and C = —2.612 are two coefficients of the expansion 



oo —nx 



E 

71=1 



n/ 



5/2 



Ax 3/2 + S + Cx + i5x 2 + ... 



(36) 



Now, using eq . (|30"D , we have the standard result 



dC- 



v 



dT 



dC^ 



T, + 



dT 



3gC 2 3.66g 



T r A 2 



T„ 



(37) 



The high temperature expansion could be obtained by using eq.([l7|), like in 



ref.[|15|- Here we shall use eq. ([18 ) , because for the aim of the present paper it 
is more convenient. Using the properties of C{s;Q±), which we have discussed in 
section 2, we see that the integrand function in eq. (|l8|) 



C R ( s + l)T(s)[C(s;Q + ) + C(s;Q-)}(3- 



-(•+1) 



(38) 



has simple poles at s = 3, 1, 0, —3, —5, —7, . . . and a double pole at s = — 1. In- 
tegrating this function on a closed path containing all the poles, we get the high 
temperature expansion, valid for T > T c (here 7 is the Euler-Mascheroni constant) 



£(/3,/i) 



■V 



71 



45/3 4 12/? 2 



a 2 - lii 2 ) + 
4 



(a 2 - /i 2 ) 3 / 2 
6tt/? 



+ £l(3<* 2 -. 2 ) + ^b(log!? + 7-i) 



2 I T7 2 ' ' 167T 2 



(39) 



■2E 



C«(2n + l)6_ (2n+1) /^ 



„ , (-27r)»(2n + 1) 

where we have used the formula 



Ci(-2n) 



T(2n+l)Ci;(2n+l) 
f-27r) n 



(40) 



Deriving eq.(|39[) with respect to fi and using eq. (|10|) we get the expansion for 
the charge density 



fiT 2 fjT{a 2 - fi 2 ) 1 / 2 fi{3a 2 -2fi 2 ) 



2tt 

T/ 2^ 



127T 2 



(41) 



V 



71=1 



;-27r) n 



For // = a, the leading term of this expression gives again the result (24). 

From eq.([|l|), by a strightforward computation and taking only the leading terms 
into account, one gets 
12tt 2 



/At; 



9« 



Td T8 (T,fx)\ 



T=T? 



-2q 



(42) 



and finally, from eqs.(|30D and (J31p 



dC- 



v 



dT 



dC^ 



dT 



32S ' 2 (43) 



9a 



in agreement with the result of ref . || . 

6 Conclusion 

We have shown that both massless and massive scalar fields on the Lobachevsky 
space H 3 exhibites Bose-Einstein condensation at a critical temperature depending 
on the curvature of the space: the higher is the curvature radius the higher is the 
critical temperature. The treatment is not intended to be complete in any sense but 
we were able to display a curvature effect on thermodynamic quantities at the most 
elementary level. In particular, due to curvature, massless charged bosons have a 
finite T c in contrast to the flat space result. The difference can be traced back to the 
existence of a gap in the spectrum of the Laplace operator on H 3 . We also pointed 
out that the infinite volume limit is under good control only for Dirichlet boundary 
conditions, for which the smallest eigenvalue has its limit value precisely at the gap 
of the continuous spectrum. 

The manifold R x H 3 with the product metric is not even a solution of Einstein 
equations nor it has constant sectional curvature. In fact, the unique simply con- 
nected static solution of Einstein equations with hyperbolic spatial part is anti-de 
Sitter space-time for which there is a non trivial goo component of the metric ten- 
sor representing a constant gravitational field. Then the equilibrium temperature is 
position dependent and a precise discussion requires the elucidation of the notion of 
local equilibrium, as has been stressed in ref.[^5|]. The same remark can be applied 
to any static gravitational field. More generally, one can ask whether there are sta- 
tionary geometries, i.e. a space-time with a Killing vector field with non-vanishing 
vorticity, such that the spatial 3-geometry is locally hyperbolic. This would be the 
gravitational counterpart of an ideal Bose gas in a rotating vessel, for which there 



are critical velocities at which correspond the appearance of quantized vortices [p7|. 
However, it has to be noted that a Robert son- Walker solution of Einstein equa- 
tions, together with the assumption that the material content of the universe ex- 



pands adiabatically, can represent a manifold of the form we have considered. The 
problem we have studied then can find physical applications in the standard model 
of the universe. 
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